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Abstract 

The vector potential for the Ohmic heating coil system of a tokamak is ob- 
tained in semi-analytical form. Comparison is made to the potential of a simple, 
finite solenoid. In the quasi-static limit, the time rate of change of the potential 
determines the induced electromotive force through the Maxwell-Lodge effect. 
Discussion of the gauge constraint is included. 



1 Introduction 

The physics of the solenoid [T] has recently enjoyed a renewed interest among investiga- 
tors [21 131 HI |5] for what it can say about the relationship between, or even the relative 
reality of, the vector potential and the electromagnetic fields. While the Aharonov- 
Bohm effect [6j [TJ, [8] has received most of the attention, its classical analogue, the 
Maxwell- Lodge effect [91 [10] suggests that electromotive induction is best described in 
terms of the potential. 

With an aim for developing a model for the electric field applied to toroidal plasma 
confinement devices, we restrict the scope of this paper to the quasi-static case appro- 
priate for low-frequency phenomena. The interested reader is directed to the paper by 
McDonald [3] and references therein for a treatment in terms of retarded potentials. 

We begin by deriving the vector potential for a simple solenoid of finite length, 
accounting for every current which one can see. The model for the circulating potential 
is built from that obtained for a single current loop. Then the Ohmic heating coil 
system of a tokamak is modelled by a circulating potential including contributions 
from the outer current loops. The use of field-cancelling outer loops can produce a 
significant region with no magnetic field present which still experiences electromotive 
induction. We conclude with a discussion of the relative merits of the potential and 
field formulations. 
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Figure 1: Current path for the simple solenoid fed from infinity with horizontal leads 
at Y = along X. 



2 Simple Solenoid 

The full expression for the vector potential of a simple solenoid turns out to be not 
so simple after all! Once one accounts for the actual path of the current, Figure [fl 
even with some continuum approximations, the result contains terms which break the 
axial symmetry (for a single feed wire). We work in (Z,X,Y) Cartesian and (Z,R,<p) 
cylindrical coordinates, where dx = d/dX, with the axis of the solenoid along Z. The 
current Iq is fed from a source at infinity along the Z-axis, and the solenoid's spatial 
extent has a height of 2h and a width of 2w with = N/h turns per unit length. 
The coils are fed from the Z-axis by leads at Y = along X such that Ix = lo- 
Along the coils, the current vector has two components, J CO ii = ZIz + such that 
I(l>/ Iz = 2nwNh and = 1 + (I(f>/Iz) 2 - For Nh 1, we make the continuum 

approximation J coi i — > K co n = I co n/2iiw, representing the line current of the conductor 
by a surface current along the cylinder with both polar and axial contributions. 

We note what appears to be the slightest of errors in the geometry used by Rousseaux, 
et al [10], in that the current path along the coil does not constitute closed loops at a 
fixed inclination to the vertical, which would break the symmetry by the polar angle 
of the inclination and does not allow for the net flow of charge along the solenoid. The 
vector potential of a tilted loop should pick up an additional rotation compared to their 
expression, as follows. Let tan#' = l/Nh2nR' define the pitch angle, as in Figure EJ A 
circular loop with radius R' at angle 6' to the horizontal midplane has |X| < R'cosO', 
\Y\ < R', and \Z\ < R'sinO'. The unit vector normal to its area has components 
(— sin 9', 0, cos 6') in (X, Y, Z) thus breaks the polar symmetry ^ in (R, <f>, Z). Us- 
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Figure 2: Coordinates for a tilted loop normal to Z . 



ing C a = cos a and similarly for sin and (X', Y', Z') for vector components in Cartesian 
coordinates relative to the tilted loop, the spherical components of the potential at a 
point P are 
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relative to a vector on the right giving the potential Ai oop . On the positive X-axis, 



and A = 4> A 
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= (jjAp, but along the positive F-axis (ft 
-XA^Cqi — ZA^Sqi = (f>A^/Co' — ZA^/So' 



ft = tt/2 and A 
dA^Se' + 4>A^Cqi. 

The expression for the magnetic field of an infinitely long, infinitely thin wire is 
well known, -B wirc = 4>hqIq/2kR. We agree with the position that the symmetry of the 
potential should match the symmetry of the constituent source, and for a static current 
density V ■ J = we use the Coloumb gauge constraint V ■ A = 0. Accordingly [TT] . 
Ampere's law becomes V 2 A = —fioJ, with solution (writing [i! Q = /i /47r) 



Ax(r) 




dr 
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for sources of finite extent vanishing to infinity, where the subscript X reminds us that 
only Cartesian components may be extracted from the integral directly, and similarly 
for surface and line currents K = 5 a J and J = 5\K. Sources which extend to infinity 
require various "tricks of the trade" for the determination of their potential; for the 
infinite wire, inspection reveals that A wire = —Zfi' Q I \nR 2 satisfies V x A wirc = B wire . 
For our current path there is a missing segment of current between points P and Q in 
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X (m) X (m) 

Figure 3: Vector potential (left column) and field (right column) produced by a simple 
solenoid with h = w = lm and iVj, = 5 turns per unit length carrying IMA of current. 
The horizontal contributions (a) and (b) are evaluated at Y = while the vertical 
and solenoidal contributions have polar symmetry. Note that the contours are linearly 
spaced for A and logarithmically spaced (base 10) for B. 



the figure, whose potential is written 

Ai g = Z/i' /o j H dZ'[(Z - Z'f + R 2 }- 1/2 , (4) 

h-Z+[{Z-hf + R^ \ 
—h — Z+ [(Z+ h) 2 + R 2 ] 1 / 2 J ' 1 ' 

and one may recover the infinite wire's expression^ by doubling the semi-infinite integral 
for the potential in the plane Z = given by taking the limit h — > oo and dropping the 

cf J. B. Tatum, Physics Topics: Electricity and Magnetism at 

http://www.astro.uvic.ca/~tatum/elmag.html 
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infinite constant's contribution, 

A™ g = Zfi' Q I 21n[z' + {Z' 2 + R 2 ) 1 / 2 

-> -Zfi' I lnR 2 = A wire . (7) 
For the horizontal feeds, one has the axial symmetry breaking terms 

A horiz = X// I x (8) 



X + [(X - wf + Y 2 + (Z- Z') 2 } 1 ' 2 



-X + [X 2 + Y 2 + {Z- Z') 2 } 1 / 2 



for Z 7 = ±h, evaluated for the plane along the X axis in Figure [3] (a) with its field 
in (b), where the potential changes sign at the horizontal midplane and the vertical 
component of the field vanishes at Y — 0. 

Next we tackle the vertical current flowing along the coils at the coil radius R' = w, 
which is equivalent to a surface current flowing along a finite, hollow pipe Kp ipe . It 
produces a potential with axial symmetry, 



A p£ = Mdb / ^' R ' / ( 9 ) 



ift.ii' _ 7 / iz f n . r dz y 

for \r — r'\ 2 = (Z — Z') 2 + R 2 + R' 2 — 2RR' cos((f) — </)'), which we wish to compute along 
the X-axis where = 0. For either order of integration, while the first is tractable, the 
second resists being put into closed form and requires a numerical evaluation. For the 
infinite pipe h — > oo, its integral 

Agf = -Z^J 7T d<P>\n(R 2 + R' 2 -2RR'cos^ , (10) 

= -Zn' — [ d0'ln(a-26cos0') (11) 

may be parametrized by a = R 2 + R' 2 and b = RR', where a > 2b with equality 
at the conductor. A direct evaluation thereof requires attention to the logarithmic 
branch which previously has led us into difficulties. Avoiding some complex analysis, 
we expand the integrand using c = 2b/ a, 

In (a - 26 cos 0') = In a + ln(l - ccos 0') , (12) 
= lna-£r=i( ccos A, (13) 

and integrate under the sum for non-vanishing even powers, 

/Vln(l-ccos0') =-EZX^/-^p^, (14) 

= -Er=ix/- d ^ c °s0') fc , (15) 
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Figure 4: Vector potential (solid), magnetic field (dashed), and current density (dash- 
dot) for an infinite pipe at increasing order of evaluation: (a) k max = 10, (b) k max = 100, 
(c) k max = 1000, and (d) k max = 10, 000. The equivalent current on axis equals 4n. 



where the beta function f3(x, y) = T(x)T(y)/T(x + y) may be expressed 
of the gamma function, to obtain the result 
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in terms 
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With normalization fi = 1 and Iz = 47r, inspection of the potential as the count of 
terms k max increases, Figure HI reveals that one may write 
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where H(x) is the Heaviside step function. Its curl -B°° p ^ = —cpdnAz gives the field 
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where the second term vanishes on account of the Dirac function S(x) = \a\5(ax), and 
its Laplacian is 

2 „,„,„, _ 2 
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whence the surface current equals K~z = 2/R' recovering 2irR'K~z = Arc = Iz- For the 
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finite pipe, one may verify that 



A h p £ = £ NtB!A«* = £ AZ'A^ (22) 



jive equivalent evaluations in the limit A<p', AZ' — >• 0, where 



'27Ti2- 



h - Z + [(Z - h) 2 + a - 26 cos 0'] 1 / 2 



In 



-/i - Z + [(Z + /i) 2 + a - 26 cos 0'] 1 / 2 



and A rin ' g will be given later. When adding up these contributions for the solenoid, 
A vert = A wire — A seg + Apipe in Figure [3] (c), one finds the relative magnitude A pipe /A seg 
is reduced by both the pitch angle Nh and the circumference 2ttR'. The corresponding 
field in (d) just barely shows the influence of A pipc beyond the radius w even for a 
contrivedly low turns per length of Nh = 5. 

Finally comes the potential arising from the solenoidal current, evaluated at <fi = 
with contributions from the Y component of the coil current 



'27T 



-h 



for \r — r'\ 2 = a — 2b cos </>', where now a = (Z — Z') 2 + R 2 + R' 2 > 2b with equality at the 
location of the loop. The angular integral evaluated for a single current loop at Z 1 may 
be expressed in terms of complete elliptic integrals with parameter k 2 = —46 /(a — 2b) 

as 



A*£ oc f ,, (25) 

> /_ 7r y/a — 2b cos <p 

= 7— L== [aK(k) - (a - 2b)E(k)} , (26) 
by a — 2b 

whereupon writing k — > ik and k' = v/I — (iA;) 2 produces a real modulus fe" = k/k', 
and the elliptic integrals transform [13] as X(i&) = K(k")/k' and E(ik) = E(k")k', 
yielding 

Aoof « 7-7=^ - (« + 2&)£(**)] , (27) 

6ya + 26 



2Va + 26 



a + 26 
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(28) 



for fc" = ^/Ab/{a + 2b), which compares well with other derivations [U E] following 
examples by Jackson [15] and Smythe [16], noting that the current around the in- 
finitesimal loop differs from the current around one of the solenoidal loops by a factor 
I'i = Sz'K^ = 5z>I ( f > /27i-R' . The remaining integral over Z 1 has so far eluded cap- 
ture in closed form; however, its numerical evaluation on a rectangular grid is simply 
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Central coil Outer coils 



Figure 5: Coil locations of a typical Ohmic heating coil assembly. Two sets of coils in a 
mirrored configuration are connected in parallel to the same current source. The outer 
coils serve to cancel the field within the confinement region. 



accomplished by adding up shifted copies of the loop potential with the appropriate 
normalization, A sol = ^2 Z , AZ'A lo ^ , for inter-loop spacing A Z i. The summation re- 
sults in the potential and field shown in Figure E] (e) and (f), where the sign of B z 
changes from the inner to outer region and Br vanishes on the horizontal midplane. 
By a similar calculation, one obtains the potential for a ring of current appearing in 
Equation (l22j) . 



ring 



-.K{k") 



'2ir^+2b~ v ~ ' ' (29) 
The final expression for the vector potential of a simple solenoid is then written as 

-^simple -^horiz -^vcrt "sol- 



3 Ohmic Heating Coil 

The Ohmic heating coil system of a typical tokamak [T7] consists of a central solenoid 
connected in series (ie with the same current) to outer loops designed to exclude the 
magnetic field from the confinement region, shown schematically in Figure |5j The 
central solenoid often consists of a doubled coil fed by a bi-directional cable, thus 
the vertical and horizontal feed currents effectively cancel, leaving only the solenoidal 
contribution. Details of the connections to and between the outer loops have not yet 
been provided, so we must here neglect their effect. We model the outer loops by two 
pairs of coils with vertical symmetry located at R' = 2.5w and 3w and Z' = ±h and 
±h/2 with Ni = 8 and jV 2 = 11 turns, respectively, and neglect their spatial extent. 
The number of turns has been selected to minimize the integral of the magnetic field's 
magnitude over the confinement region, which we approximate with a circular vessel of 
radius w centered at R = 2w. The magnitude of the current swing can be estimated as 
±110kA per second [17], with a pulse duration exceeding several seconds. The electric 
field distribution on the order of mV/m may be read directly from the vector potential 
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Figure 6: Potential (a) and field (b)-(d) produced by the solenoidal component of 
a typical Ohmic heating coil system carrying 11 Ok A of current with h = w = lm, 
Nh = 50, and two pair of outer coils. The white circle represents the confinement 
vessel, at whose center the field vanishes but the potential does not. Note that the 
contours are linearly spaced for A and logarithmically spaced (base 10) for B. 

in the quasi-static approximation i" — > d t I as E = —d t A when V$ = 0. 

The presence of the outer coils serves to flatten the potential as shown in Figure[6](a), 
thus reducing the magnetic field within the confinement region, as seen in (b). Both 
the vertical (c) and horizontal (d) field components change sign within the vacuum 
vessel, and for a sufficiently complicated outer coil configuration, a field-free region of 
significant extent can be achieved. The exclusion of the field in no way affects the 
electromotive induction present at the center of the vessel resulting from the time rate 
of change of the vector potential <9tA h c , which does not vanish within the confinement 
region. The situation begs the question, how does a charged particle know that it 
should move when at the location of vanishing magnetic field? We feel the answer lies 
in ascribing to the classical potential a sense of reality which surpasses that of the field 
formulation. 

We close this section by noting that the contours of do not give directly what 
would be called the flux surfaces for this configuration. For KV ■ Bzr = with = 0, 
the field B Z r may be written as the contours of a flux function {—Or, dz)ip = 4> x V?/> = 
>kip such that = RB Z r = -RV x A^, thus W(RA^ + ip) = and the flux function 
is given (up to an unphysical constant) by — RA^. The central solenoid's contribution 
is shown in Figure [7] (a) with the field's angle from the vertical computed directly from 
tan 6* = Bz/Bn. Compared to the flux function of the simple solenoid, the Ohmic 
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Figure 7: Contours (solid) of the flux function ifrip = RBzr for —ip = RA^ give the 
field lines for the central solenoid (a) and the Ohmic heating coil system (b). Also 
shown are contours (dotted) of the field line's angle evenly spaced in units of 7r/4. 

heating coil system (b) has a much more complicated topography induced by the outer 
coil assembly. 

4 Discussion and Conclusion 

The argument over which is more fundamental, the potential or the field, has yet to be 
decided [HI QU EH EO [221 US [H] . From the hydrodynamical analogue [25], one sees 
that the choice of (Coulomb)Lorenz gauge reflects a statement on the (in) compressibility 
of the electromagnetic potential/flow and on whether the speed of propagation of dis- 
turbances is (in)finite. Accordingly, in situations when V • J = 0, one should take 
the Coulomb gauge V ■ A = 0, and when a varying space charge density may develop 
dtp e 7^ 0, then the Lorenz gauge d^A^ = is appropriate. From the geometric state- 
ment [261 123 (2H1 [29] of classical electromagnetism in terms of dual forms d*dv4 = J, an 
intimate relationship is apparent between the potential and the source. We feel that 
this relationship is beyond cosmetic and speaks to the fundamental structure of natural 
phenomena. 

In executing the solution of the Poisson equations Q, one evaluates at a point 
(Z, X, Y) the weighted contribution from all source points (Z', X', Y') which have had 
time to communicate their influence to the point in question, which for the Coulomb 
gauge with infinite propagation speed of disturbances extends to infinity. Conversely, 
the action of the current source has made an image of itself which disperses through 
the region of consideration at the speed of disturbances, which we call the potential. 
The description in terms of diverging and circulating fields adds a level of abstraction 
analogous to ascribing to the slope of a hill the same reality as the hill itself; while the 
slope of the hill informs a ball how to roll under gravity, it is the material of the hill 
which is real. Note that we are not advocating a revival of the luminiferous aether |3Uj . 
but simply a more fundamental view of the nature of the potential. Its supposed 
indeterminacy reflects the fact that a particular set of fields within a region can be 
generated by a multitude of source configurations external to or within that region; 
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however, once the source has been specified, then so has the shape of the potential, 
hence the fields. Just because we cannot see it does not mean it is not there. 

Operation of a tokamak (or similar, smaller apparatus [10]) affords the opportunity 
to explore phenomena beyond that associated with fusion and to think about issues of 
basic science, such as the nature of induction. The description of complicated current 
configurations becomes straightforward in the potential formulation. Examination of 
the vector potential applied to a tokamak indicates that it is a prime example of the 
Maxwell-Lodge effect, where one tries to understand how a charge knows it is to move 
when there is no field changing with time at its location. 

Attempts to ascertain what is happening within the device [31] would do well to 
begin by determining what is being done to the device. The true starting point of any 
calculation involving tokamak operation is the determination of the potential produced 
by sources external to the confinement region. The usual evaluation in terms of resistive 
magnetohydrodynamics utilizing the quasi-neutral approximation requires the electric 
field be determined from an equation of motion rather than Poisson's equation for the 
potential, thereby denying any predictive use of the model. From a determination of 
the applied potential, one should in principle use the Ohm's law equation (conspicuous 
by its absence beyond the kinetic term in our treatment [32] of the macroscopic field 
formulation) to predict the resulting current density. Only then can one treat both the 
fluid and the potential through a hydrodynamic field theory. The goal of experimental 
science is to predict results, not to interpret them, and for a tokamak that requires 
modelling the electric field which the investigator has applied. 
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